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Abstract 

We present a time-resolved terahertz spectroscopic study of the half-metallic ferromagnet Cr02. 
The ultrafast conductivity dynamics excited by an optical pump displays a very short (several 
picoseconds) and a very long (several hundred picoseconds) characteristic time scales. We attribute 
the former to the electron-phonon relaxation and the latter to the spin-lattice relaxation. We use 
this distinction to quantify the relative contribution of the scattering by spin fluctuations to the 
resistivity of Cr 02 ; we find that they contribute less than one half of all scattering events below 
room temperature. This contribution rises to ~ 70% as the temperature approaches Tc=390 K. 
The small effect of spin fluctuations on the resistivity is unexpected in the light of the proposed 
double-exchange nature of the electronic and magnetic properties of Cr 02 . 
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I. INTRODUCTION 


Chromium dioxide Cr02 is a half-metallic ferromagnet (Tc=390 K), in which the ma¬ 
jority spin electrons are metallic while the minority spin electrons are semiconducting, i.e., 
the Fermi level falls within a gap in the minority density of states^. The nearly 100% spin 
polarization^i^- makes Cr 02 attractive as the source of spin-polarized electrons in spintron- 
ics, while the material was also used as the magnetic recording medium. Theory predicted^ 
the magnetic moment per Cr^"*" (Sd^) ion to be 2/is, in agreement with Hund’s rules and 
experiment^. Of the two d electrons, one is localized and found about 1 eV below the Fermi 
level. The other d electron hybridizes with the oxygen p orbitals and forms a narrow itinerant 
band that crosses the Fermi level. Korotin et ah- used the term ”a self-doped double ex¬ 
change ferromagnet” to describe the material’s intertwined metallicity and ferromagnetism, 
with the mobile d electrons mediating the double exchange between the localized d spins. 
Another remarkable feature is the two-order-of-magnitude drop in resistivity between 400 K 
and 10 K (Fig. [2]) whose origin is not fully understood. In this work, we use time-resolved 
terahertz spectroscopy (TRTS) to compare the relative roles of spin fluctuation (or spin-flip) 
and other scattering processes in the resistivity of Cr02. We hnd that the spin-flip processes 
do not dominate the electron scattering in a wide range of temperatures below Tc, as many 
authors have assumed. 

Cr02 crystallizes in the tetragonal rutile structure, with lattice parameters a = b = 0.4421 
nm and c = 0.2916 nm^. The Cr atoms are octahedrally coordinated by oxygen, and edge¬ 
sharing oxygen octahedra form ribbons along the c axis, while the octahedra on adjacent 
ribbons share a corner—. The Fermi level for the majority spins belongs in a half-hlled 
band derived from the the dxy and dyz orbitals^^l^. The band gap in the minority density 
of states exceeds 2 eV, with the empty minority states lying about 1 eV higher than the 
Fermi level^^i^. The half-metallicity of Cr02 was conhrmed experimentally by point contact 
Andreev spectroscopy^, tunneling measurements^, and spin-polarized photoemission^. 

A survey of literature hnds no agreement on the origin of the temperature dependence of 
resistivity (Fig. [2]). Lewis et ah— showed that below about 200 K, the temperature depen¬ 
dence is well described by the Bloch-Gruneisen function and phonon scattering dominates. 
Above 200 K, spin-flip scattering becomes important and contributes one half of all scat¬ 
tering events near the Curie temperaturel^i^i. Barry et ah— suggested a phenomenological 
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description based on the formula p{T) = po + AT^exp(—A/T) with a gap A 80 K, 
above which the resistivity follows the dependence expected for a spin-flip scattering in 
a metallic ferromagnet. Gupta et ah— £t the low-temperature resistivity (below 40 K) with 
a p{T) = Po + AT^ dependence characteristic of spin-flip scattering if the non-rigid band 
behavior of the minority band is accounted for. Watts et al. proposed a two-band picture 
for electronic conduction based on a magnetotransport studyl^, although other magnetore¬ 
sistance studies have not reached the same conclusiortl^'i^. Several authors found that a 
dependence also describes well the resistivity data in a broad temperature range and 
attributed this to electron-electron scattering^ii^. 

The contradictory scenarios proposed for electron conduction in Cr02 perhaps reflect the 
reality that all three scattering processes - electron-electron, electron-phonon, and spin-flip - 
play a role. Our TRTS study is motivated by the possibility of separating the different scat¬ 
tering contributions based on the different time scales for the coupling of electrons, lattice, 
and spins to the optical pump excitation. When a metal is excited by the optical pump, 
the absorbed photons deposit their energy into the electronic system. Within a picosecond, 
the relaxation of this energy establishes a thermal electron and phonon distribution at an 
elevated temperature^^'— . The subsequent thermalization of spins happens much slower, on 
the scale from tens to hundreds of picoseconds^ii^. This vastly slower spin thermalization 
allows us to distinguish the contribution of the spin-flip scattering to resistivity from the 
contributions of the electron-electron and electron-phonon scattering. 

II. EXPERIMENTAL DETAILS AND RESULTS 

For this study, we used an epitaxial 100 nm Cr02 thin him grown on a (lOO)-oriented 
0.5 mm Ti02 substrate using chemical vapor deposition with CrOa as precursor—. TRTS 
and terahertz time-domain spectroscopy (THz TDS) measurements were performed using 
a home-built spectrometer based on an amplihed Tiisapphire laser with 1 kHz repetition 
rate^^. The THz wave was polarized along the crystalline b axis of Cr02. The optical 
pump pulses with 800 nm wavelength (1.55 eV photon energy) and 0.2 mJ/cm^ huence were 
polarized along the c axis. The THz probe spot diameter was 2 mm, while the diameter of 
the optical pump spot was 10 mm. Temperature control in the 9-400 K range was provided 
by a He how or a closed cycle cryostat. 
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Time-resolved magneto-optical Kerr effect measurements (MOKE) were carried out at 
room temperature in the polar MOKE configuration. The pump and probe wavelength was 
800 nm. The static applied magnetic held was normal to the him and measured to be 3300 
Gauss. It was supplied by a stack of permanent magnets. 



Frequency (THz) Frequency (THz) 


FIG. 1: (Color online) (a) Red line: THz pulse transmitted by the Cr02 film in equilibrium state 
(no pump excitation) at 9 K. Blue line: the transmitted THz pulse 200 ps after the optical pump 
at 9 K. (b),(c) Real and imaginary parts of the THz conductivity at 9 and 70 K in equilibrium 
state. Symbols: measured conductivity. Solid lines: the Drude model fit. (d) THz conductivity 
at 150 K and 330 K. Only the real part is shown for the 330 K data. Solid lines show the Drude 
model fit for the 150 K data. 


THz TDS measures the electric held of a THz pulse in time domain. Figure [T](a) shows the 
THz pulse after passing through the Cr02 him at 9 K. To extract the THz conductivity of the 
him, a bare Ti02 substrate was used as reference^^. The sample and reference measurements 
are Fourier-transformed to obtain the frequency domain spectra Sgami^^) and S're/(n;) and 
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compute the amplitude transmission coefficient t{u) = Ssami,^)/Srefi,^)- We compute the 
THz optical conductivity from2ii2^ 


~ ha + l [ .u:{ds - dr){n‘i - 1) 

= V- ry _/ ^ J 6Xp ( I 


( 1 ) 


hs + ZQa{oj)d 

where 77,3 is the THz refractive index of Ti 02 , d is the film thickness, Zq = 377 H is the 
free space impedance, and {ds — dr) is the difference in thickness between the film and the 
bare reference substrates. The frequency dependence of the conductivity is well described 
by the Drude model cr^u) = o'o/(l + iuj/'y), where Uo is the dc conductivity and 7 is the 
electron scattering rate, the parameters whose temperature dependence is determined by 
least-square fits (Fig. [T]^b-d)). The real conductivity becomes frequency-independent above 
150 K in our THz frequency window and the scattering rate 7 is not reliably measured. 
The temperature dependences of Fig. |2] agree well with the published transport and optical 
conduct ivity 2 ^ studies: cxo undergoes a two-order-of-magnitude change, while 7 also drops 
precipitously to ~ 0.5 THz at 9 K. A similar ’’collapse of the scattering rate” was found by 
Singley et al. and is responsible for the low residual resist ivity2^. 
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FIG. 2: (Color online) The temperature dependence of the dc resistivity and the scattering rate 7 
determined by THz TDS. Drude fit error bars are smaller than the size of the symbols. 

The THz pulse lasts only several picoseconds (Fig. [11(a)), which allows the measurement 
of conductivity changes with picosecond time resolution. Figure [T](a) shows the transmitted 
THz pulses before and 200 ps after the him is excited by a sub-lOO-fs optical pump pulse. 
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The THz probe pulse transmitted before the pump measures the equilibrium conductivity; 
the probe that passes after the pump measures the conductivity in a non-equilibrium state. 
At 9 K, the non-equilibrium THz pulse displays a higher electric field amplitude and a 
different phase relative to the equilibrium pulse (Fig.[T]^a)). Above 90 K, the phase difference 
between equilibrium and non-equilibrium THz probes becomes immeasurable, while the THz 
amplitude remains higher for the non-equilibrium probe. The higher amplitude and the 
different phase of the transmitted THz probe result from a lower conductivity and a higher 
scattering rate in the non-equilibrium state. These pump-induced conductivity changes 
indicate an elevated instantaneous temperature in the evolving non-equilibrium state. The 
pump-induced changes are consistent with the hndings of Fig.[2]that show lower conductivity 
at higher temperature. Thus, the effect of the optical pump is a very fast, picosecond-scale 
heating of the Cr 02 him. 

By varying the arrival time of the THz probe relative to the optical pump, we record 
the evolution of pump-induced conductivity changes. For simplicity, we only measure the 
change in the peak transmitted THz held instead of recording the full THz pulse. Figure [3] 
shows the measured relative change AE{t)/Ef) as a function of time delay between the pump 
and probe; Eq is the peak THz held in the absence of the pump excitation. 

Since conductivity a{oj) is (almost) independent of frequency oj above 150 K, we can take 
cr^u) ~ (To and relate the change AE{t)/EQ to the frequency-independent pump-induced 
change A(T as 

AE{t) ^ -ZodAa 

Eq l + h^ + Zodao 

On the right hand side of Eq. ([2]), only Aa contains the ehect of the optical pump. All 
other quantities characterize the equilibrium state. Thus, the time-evolution of AE{t)/EQ 
rehects the time-evolution of Aa. A higher transmitted THz held (positive AE) indicates a 
drop in conductivity (negative A(t). 

Figure [3](a) shows two vastly diherent time scales in the conductivity response to the 
optical pump. A fast step-like rise in AE{t) is followed by a much slower evolution, as 
AE{t) reaches a broad maximum (near 400 ps at 300 K in Fig. [3](a)) and then begins a slow 
recovery of its equilibrium value. Below 250 K, the initial rise in AE{t) is followed by a fast 
shallow drop (Fig. [3](b)) before the broad maximum and the recovery of equilibrium. The 
broad maximum in AE{t) is found at all temperatures down to 70 K but becomes a lot less 
pronounced below 150 K (Fig. [3](a)). 
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FIG. 3: (Color online) (a) TRTS spectra - the time-resolved measurement of the ultrafast change 
in the peak transmitted THz electric field at various temperatures, (b) TRTS spectra similar to 
(a), but zoomed in on the first 30-60 ps of the photoinduced response and recorded with a shorter 
time step. The spectra were shifted horizontally for clarity. 


III. DISCUSSION OF THE RESULTS 

The dynamics of AE{t) reflects the relaxation of the pnmp excitation energy. The 1.55- 
eV pnmp photons are absorbed by the transitions in the majority channel, as the gap 
in the minority channel exceeds 2 eV. The optical pnmp creates a highly excited non- 
thermal electron popnlation. The evolntion of this excited state is nsnally described in 
terms of a fast (~ 100 fs) electron thermalization at an elevated temperatnre, which is 
followed by the electron-phonon relaxation and the eqnilibration of electronic and phonon 
temperatnresl^‘22i2ii2^i2^. Zhang et al. fonnd that at 300 K in Cr02, the electron and phonon 
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temperatures reach an equilibrium in about 2-3 ps^, which is consistent with our data. The 
fast shallow drop in AE{t) at T < 250 K corresponds to a slight recovery of conductivity as 
energy is transferred from electrons to phonons (Fig. [3]^b)). The absence of this conductivity 
recovery feature at high temperature indicates that the phonon scattering and the phonon 
temperature gain relative importance in the conductivity dynamics. With or without the 
slight recovery, we interpret the initial (< 5 ps) dynamics in AE{t) as the electron-phonon 
relaxation, after which elevated and equal electron and phonon temperatures are established, 
leading to higher electron-electron and electron-phonon scattering rates. 

Why does the conductivity Aa{t) continue to drop (the resistivity Ap{t) continue to 
rise) after the initial dynamics? Another process that contributes to resistivity is the spin- 
flip scattering. The evidence for the spin temperature evolution in Cr02 after the optical 
pump is provided by the magneto-optical Kerr effect^^^^ (MOKE), which refers to a change 
in the polarization state of reflected light and is proportional to the material’s magnetiza¬ 
tion. In time-resolved MOKE (TRMOKE), the pump-induced change in magnetization is 
recorded^^— . In Cr02, a slow demagnetization over hundreds of picoseconds follows the 
optical pump excitation, as the spin temperature rises due to the spin-lattice couphng2Ii2^. 
Figure 111(a) shows a room-temperature TRMOKE measurement in which we recorded the 
pump-induced polarization rotation of an optical 1.55-eV probe pulse; we observe a fast 
initial jump and a much slower increase over 1 ns. The comparison with published data^^*^^ 
shows that the slow TRMOKE dynamics reflects the demagnetization as the spin temper¬ 
ature equilibrates with the electron and phonon temperature via the spin-lattice coupling. 
The room-temperature spin-lattice relaxation time was measured by Zhang et al.— to be 
~ 400 ps. Thus, the broad maximum in AE{t) results from the rise in the spin temperature 
and the corresponding increase in the spin-flip scattering. 

To quantify the effect of the spin fluctuations on resistivity, we turn to the two- 
temperature model of the coupled electron-lattice and spin dynamics^*^. This model is 
applicable for the long time-scale (f > 5 ps) dynamics, when the electron and phonon tem¬ 
peratures can be taken as equal and described by a single electron-phonon temperature Tep. 
The spin system is described by the spin temperature Tg. The spin and electron-phonon 
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FIG. 4: (Color online) (a) TRTS and TRMOKE spectra at room temperature. Solid lines are fits 
to the two-temperature model, (b) TRTS spectra at 150 K and 380 K. The solid lines are fits to 
the two-temperature model. 


dynamics are described by a pair of differential eqnations 


dT, 


ep 


•'ep 


dt 


dTs 

' dt 


-t/T)]Tep-Tb) , 

( 3 ) 

= -9 (Ts - Tep ), 

( 4 ) 


where Cgp and Cg are electron-phonon and spin specific heats and g is the spin-lattice conpling 
constant^. The specihc heats c^p and Cg are taken as temperatnre indepedent nnder the 
assnmption of a small pnmp-indnced temperatnre change in the electron-phonon and spin 
systems. The last term in (|3]) describes the cooling of the electron-phonon system by the 
diffnsion of energy into the snbstrate, whose temperatnre Tf, is taken as constant and eqnal 
to the eqnilibrinm temperatnre of the measnrement. The cooling is proportional to the 
temperatnre difference (Tgp — Tt) and is parametrized by a constant /3. The exponential 
that mnltiplies T^p in the last term of Eq. ([3]) acconnts for the gradnal ’’tnrning on” of the 
cooling as the energy deposited within the optical absorption depth of the pnmp wavelength 
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diffuses through the him thickness to the substrate side of the him. We estimate the optical 
penetration depth to be 17 nm using the optical constants measured by Stewart et al.— 
The electron-phonon specihc heat c^p consists of the electronic and phonon specihc heats. 
We compute the electronic contribution as^ Cg = 7 T with 7 = 7 mJ K“^mole“^. We 
compute the acoustic phonon specihc heat in the Debye model with 6^ = 593 The 
optical phonon specihc heat is computed in the Einstein model using the frequency Uopt = 450 
cm“^ to represent the branches of the optical phonon spectruirt^S. For the spin specihc 
heat, we use the mean-held value^I Cs(T) = —where a = 3SRTc/[2{S -|- l)M(O)^], 
S = 2 for Cr 02 , and M(0) is the saturation magnetization at low temperature. We use the 
magnetometry measurements of Li et al.— to determine M(0)^ and ^r. The computed 
electron-phonon and spin specihc heats are shown in Fig. [5](a). 



FIG. 5: (Color online) (a) The computed electron-phonon and spin specific heat, (b) Temperature 
dependence of the parameter a, which quantifies the relative contribution of the spin temperature 
change to the photoinduced resistivity change. The vertical line indicates the Curie temperature 
Tc = 390 K. 

Next, we use the two-temperature model (Eqs. fl3]IT|) ) to ht the experimental TRTS data. 
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We compute the time-resolved ultrafast change in resistivity as 

f(t) = p„(r) + [(1 - a)Ar,,(«) + aAr.(«)l, (5) 

where poiT) is the equilibrium temperature-dependent resistivity, dpo(T)/dT is the slope of 
Po{T), and both are taken from Fig. [2l The parameter a describes the relative importance 
of the evolving photoinduced change in electron-phonon and spin temperatures, ATep(t) and 
A Ts(t), in the determination of the photoinduced resistivity p{t). We set the initial spin 
temperature as equal to the equilibrium temperature of the measurement, Ts{0) = T^. At 
300 K, we have both TRMOKE and TRTS data and we £t both of them simultaneously 
by assuming that the TRMOKE angle is proportional to ATs{t). Fitting both data sets 
with the same model parameters allows us to determine the parameters g and a simulta¬ 
neously. We find the value of the spin-lattice coupling constant g = 0.011 J/(mole K ps), 
which compares well to g = 0.018 J/(mole K ps) deduced from the data of Zhang et al.— 
The fractional photoinduced magnetization change was estimated to be 0.6% at 1000 ps 
time delay (Fig. 111(a)). To £t the TRTS data at the other temperatures, we consider the 
spin lattice coupling g as independent of temperature^^. Figure [5](b) shows the obtained 
temperature dependence of the htting parameter a. We tested how stable our htting pro¬ 
cedure was under the variation of the htting parameters [3 and r that describe the cooling 
of the electron-phonon system by energy diffusion into the substrate via the empirical last 
term in Eq. ([3]). We varied r between 10 ps and 10000 ps and obtained similar quality 
hts to the experimental data. While the htting parameter (3 needed to be adjusted by a 
signihcant amount to accommodate the large range of r, the htted values of a changed only 
very slightly, as rehected by its standard deviation (the error bars) reported in Fig. [5](b). 
Figure [6] shows the computed evolution of the temperature diherence (Tep(t) — Ts(t)) and 
indicates that the photoinduced instantaneous change in electron and spin temperatures is 
small at all but the lowest (T < 77 K) temperatures in our measurement. 

Boltzmann transport theory describes the electric current by summing up the electron 
velocities over the occupied quasiparticle states. The resistivity is caused by quasiparticle 
scattering between crystal momentum states, and is proportional to the electron scatter¬ 
ing rate 7 : poiT) = '-y/{eoUp), with cq being the permittivity of free space and Up being 
the plasma frequency. When the various scattering mechanisms are independent, they are 
combined using Matthiessen’s rule as 7 (T) = 7 ep(T) + IsiT), where 'yepiT) includes the 
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FIG. 6: (Color online) Time evolntion of the difference in the instantaneons electron-lattice and spin 
temperature computed in the two-temperature model. Colors correspond to different equilibrium 
temperatures. 

electron-phonon and electron-electron scattering and 7 s (7") describes the spin-flip scatter¬ 
ing. Since the plasma frequency Up in Cr 02 is independent of temperature^ below Tc, 
the photoinduced change in resistivity happens because the elevated instantaneous electron- 
phonon and spin temperatures modify the respective scattering rates, 7 ep(Tep) and ■js(Ts). 
The parameter a introduced in Eq. ([5]) quantihes the relative contribution of the spin-flip 
scattering to all scattering events. According to Fig. |5](b), the spin-flip processes do not 
dominate the electron scattering in the wide temperature range 70-300 K where the most of 
the temperature-induced change in resistivity occurs below Tq (Fig. |2]). The exception are 
the temperatures above and near Tc, where a reaches about 70% at 380 K. 

It is instructive to compare our Endings with a study by Averitt et ah— of the double 
exchange manganites Lao. 7 Ca(Sr)o. 3 Mn 03 . The manganites exhibit an ultrafast conductivity 
response with two distinct time scales, with the short ~ 2 ps scale attributed to electron- 
phonon relaxation, and the longer scale of tens of ps attributed to spin-lattice relaxatiort^S. 
The phonon-induced conductivity change Ap(ATep) dominates at low temperature (T < 
O.STc), while the spin-fluctuation-induced Ap{ATs) dominates close to Tc- This behavior 
is similar to our Endings in Cr 02 . Spin fluctuations are greatly enhanced near Tc, and 
the relative importance of the spin-flip scattering rises in both materials. However, there 
is also a significant difference between Cr 02 and the manganites: in the manganites, the 
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resistivity is highly sensitive to temperature near Tc and even exhibits a metal-insulator 
transition driven by the double-exchange physics^^. In Cr02, the resistivity is featureless 
near the ferromagnetic phase transition^^*^ (Fig. H]), despite the relative enhancement of 
the spin-flip scattering close to Tc- Our results point to a significant disconnect between 
charge transport and magnetic order, which conflicts with the double exchange scenario 
of magnetism in Cr02. Other evidence for such conflict is provided by the spectroscopic 
ellipsometry study of Cr02 films by Stewart et al.—, who found that the Drude plasma 
frequency and the effective number of carriers remain constant across the ferromagnetic 
phase transition. By contrast, in the double-exchange manganites, a significant transfer 
of the spectral weight is found from high to low energy and the Drude response grows in 
strength as the temperature is lowered below Tc^^—- 

IV. SUMMARY 

We have presented a TRTS study of the half-metallic ferromagnet Cr02, where the ultra¬ 
fast resistivity response is governed by the electron-phonon and spin-lattice relaxation. In 
the first 2-5 ps after the optical pump pulse, we observe a step-like change in the photoin- 
duced time-resolved resistivity p{t), which we ascribe to the establishment of an elevated 
electron and phonon temperature (Fig. 0]). The fast step-like feature is followed by a con¬ 
tinued slow rise in p{t) before the recovery of the equilibrium state begins. The slow rise 
in p{t) after the initial fast dynamics can be explained as the heating of the spin system 
via spin-lattice coupling. We use a two-temperature model of coupled electron-phonon and 
spin dynamics to quantify the contribution of the spin temperature change to the change in 
resistivity. As shown in Fig. [5](b), the spin fluctuations provide the dominant contribution 
to photoinduced resistivity only near Tc- Below room temperature, the electron-phonon 
and electron-electron scattering dominates. This finding should provide further guidance to 
theoretical descriptions of electronic transport in Cr02. Many of the past theoretical models 
neglected either the spin-fluctuation or the electron and phonon scattering, while our results 
show that both must taken into account. 
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